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• Abstract: We derive the allowed ranges of the spin-independent interaction cross section 

G~ for the elastic scattering of neutralinos on proton for wide ranges of parameters of the 
general Minimal Supersymmetric Standard Model. We investigate the effects of the lower 
limits on Higgs and superpartner masses from colliders, as well as the impact of constraints 
from b — > 57 and the new measurement of (g — 2) M on the upper and lower limits on g^ 1 . 
We further explore the impact of the neutralino relic density, including coannihilation, 
and of theoretical assumptions about the largest allowed values of the supersymmetric 
parameters. For /x > 0, requiring the latter to lie below 1 TeV leads to g^ > 10 -11 pb at 
m x ~ 100 GeV and g^ > 10 -8 pb at m x ~ 1 TeV. When the supersymmetric parameters 
are allowed above 1 TeV, for 440 GeV < m x < 1020 GeV we derive a parameter-independent 
lower limit of > 2 x 10~ 12 pb. (No similar lower limits can be set for fi < nor for 
1020 GeV <rn x < 2.6 TeV.) Requiring £l x h? < 0.3 implies a parameter-independent upper 
limit m x < 2.6 TeV. The new e + e~-based measurement of (g — 2)„ restricts m x < 350 GeV 
at lo" CL and m x < 515 GeV at 2g CL, and implies /i > 0. The largest allowed values of 
<7p have already become accessible to recent experimental searches. 
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1. Introduction 

The hypothesis of the lightest neutralino x, as the lightest supersymmetric particle (LSP), 
providing the dominant contribution to cold dark matter (CDM) in the Universe, has 
inspired much activity in the overlap of today's particle physics and cosmology. It is 
well-known that the relic density of the neutralinos is often comparable with the critical 
density Q g|. The expectation that the Galactic dark matter (DM) halo is mostly made 
of weakly-interacting massive particles (WIMPs) has further led to much experimental 
activity. In particular, the experiments looking for CDM WIMPs elastically scattering off 
underground targets have recently set limits on spin-independent (SI), or scalar, cross sec- 
tion of the order of 1CP 6 pb ||, ||, ||. They have also nearly ruled out the region of (m x , (Tp ) 
that has been claimed by the DAMA experiment to be consistent with an annual modula- 
tion effect || . Initial and early studies |7|, ||, ||, [R]] were followed by more recent work [|1J , 
where it was concluded that current experimental sensitivity is generally comparable with 
the ranges expected from the neutralino WIMP in the Minimal Supersymmetric Standard 
Model (MSSM). It is, however, still at least on order of magnitude, or so, above the ranges 



predicted by recent analyses of the Constrained MSSM (CMSSM) 0, [L§ . 

For comparison, cross sections for spin-dependent (SD) interactions are in the case of 
the neutralino generally some two or three orders of magnitude larger than the SI ones. On 
the other hand, at present detectors are still not sensitive enough to explore the parameter 
space of the MSSM, despite recent progress [jl7) ]. 

In light of the ongoing and planned experimental activities, it is timely to conduct a 
thorough and careful re-analysis of the predicted cross sections for SI scattering of neu- 
tralino WIMPs. Such a study is rather challenging because resulting ranges often strongly 
depend on a given SUSY model and on related theoretical assumptions. They are further 
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affected by experimental limits on SUSY, both from colliders and from indirect searches, 
as well as by cosmological input, where the relic abundance of the CDM has been mea- 
sured with better accuracy both directly and in CMBR studies [Q. Over the last few 
years and months there have been also new results for LEP lower bounds on the masses 
of the lightest Higgs and electroweakly-interacting superpartners, Tevatron lower limits 
on strongly-interacting superpartners, as well as limits on allowed SUSY contributions to 
b — > 57, and especially to the anomalous magnetic moment of the muon (g — 2) M |1S|] . The 
new result for {g — 2)^ indicates a sizable deviation from the Standard Model prediction, 
whose value is still a subject of much discussion. As we will show, when interpreted in 
terms of SUSY, the required extra contribution to {g — 2)^ plays a unique role in implying 
a stringent upper bound m x < 350 GeV (1 a CL) and m x < 512 GeV (2 a), but it does not 
affect much the allowed ranges of the SI scattering cross section a^ 1 . 

In this paper we carefully study the impact of the above constraints. In an attempt to 
minimize theoretical bias, we work here in the context of the general MSSM, which will be 
defined below, with an additional assumption of -R-parity conservation. We focus here on 
the SI cross section case. Other recent studies of the general MSSM include 
The results presented here show that the level of experimental sensitivity that has recently 
been reached ||, ||, || has now indeed allowed one to start exploring cosmologically favored 
ranges of the neutralino WIMP mass and SI cross section. However, we point out a 
number of caveats and relations and further discuss the origin of, and robustness of, the 
upper and lower limits on . In particular, for a big range of the heavy neutralino mass 
440 GeV < m x < 1020 GeV GeV we are able to derive parameter-independent lower bounds 
on <Jp. Collider lower limits on Higgs and superpartner masses plus requiring Q x h 2 < 0.3 
alone leads to a parameter-independent upper bound m x < 2.6 TeV. 



20, m m- 



2. The MSSM 

We start by reviewing relevant features of the general MSSM p3|. (We follow the conven- 



tion of [24].) By the MSSM we mean a supersymmetrized version of the Standard Model, 
with Yukawa and soft SUSY-breaking terms consistent with imparity. We neglect CP- 
violating phases in the Higgs and SUSY sectors and assume no mixings among different 
generations of squarks and sleptons, since both are probably small. In the MSSM, all the 
slepton and squark masses can be considered as a priori free parameters set at the elec- 
troweak scale. In the same way we treat the trilinear parameters A{ (i = t, b, r) of the third 
generation while neglecting the ones of the first two. The Higgs sector is determined at the 
tree level by the usual ratio of the neutral Higgs VEV's tan = vt/vf, and the mass of the 
pseudoscalar ttia- In computing full one- loop and leading two-loop radiative corrections 



to the lightest scalar Higgs we use the package FeynHiggsFast (FHF) [25]. As we will see, 
Higgs masses will play an important role in the analysis. 

The lightest neutralino x is lightest of the four mass eigenstates of the linear combi- 
nations of the bino B, the wino W3 and the two higgsinos H® and 



X = Xi = N n B + iYi 2 IU 3 ° + N 13 H$ + N 14 H°. (2.1) 
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The neutralino mass matrix A4 [£4|] is determined by the U(l)y and SU(2)l gaugino 
mass parameters M\ and M%, respectively (and we impose the relation M\ = | tan 2 6\yM2, 
which comes from assuming gaugino mass unification at GUT scale), the Higgs/higgsino 
mass parameter /i, as well as tan/3. In the region \fj,\ 3> M\, the lightest neutralino is 
mostly a bino with mass m x ~ M\. In the other extreme, it is mostly higgsino-dominated 
and m x ~ \fi\. 

For the purpose of this analysis, we take as independent parameters: tan/3, fi, M2, 
rriA, A ti b as well as the soft masses of the sleptons and of the squarks. In order to make 
our analysis manageable, we make an additional assumption that, at the electroweak scale, 
the soft mass parameters of the sleptons are all equal to some common value fhr, and 
analogously fhg for all the squarks. One normally expects certain relations among the 
physical masses of the sleptons and the squarks since they in addition receive well-defined 
D-term and F-term contributions to their mass matrices. Assuming common soft mass 
terms, at either GUT or electroweak scale, this normally leads to the sleptons being lighter 
than the squarks. Furthermore, for the sfermions of the 3rd generation it is natural to 
expect large mass splittings. However, we believe that, for our purpose, introducing just 
two separate common soft mass scales rhj and fhg, while greatly simplifying the analysis, will 
not play much role in our overall conclusions for the SI cross sections. (This is in contrast 
to often assumed full degeneracy of soft sfermion masses which leads in our opinion to an 
unnecessary limitation on the allowed SUSY parameter space.) 

What we do find important is to disentangle squark and slepton masses. This is 
mainly because experimental limits on slepton masses are significantly weaker than in the 
case of squarks. For the case of the bino-like neutralino, which is the most natural case 
for providing fl x h 2 ~ 1 |26fl , it is therefore the mass of the lightest slepton which often 
predominantly determines £l x h 2 . Assuming common soft masses for sleptons and squarks at 
the electroweak scale would therefore generally lead to overestimating Sl x h 2 by missing the 
cases where relatively light sleptons (below current squark mass bounds) would otherwise 
provide acceptable Q x h 2 . An additional effect is that of coannihilation. When slepton mass 
is only somewhat larger than that of the LSP, the neutralino relic abundance is strongly 



reduced and otherwise forbidden cases become allowed [27]. 
3. The Spin— Independent Cross Section 

For non-relativistic Majorana particles, like the neutralino WIMP, the elastic scattering 
off constituent quarks and gluons of some nucleon %X is given by an effective differential 
cross section [jj], ||, 10 1 



da da SI da SD , , 



d\q\ 2 d\q\ 2 d\q\< 

where the transferred momentum q = fj,^v depends on the velocity v of the incident WIMP, 
and f/,j\ = mAm x / (rriA + m x ) is the reduced mass of the system. The effective WIMP- 
nucleon cross sections a SI and a SD are computed by evaluating nucleonic matrix elements 
of corresponding WIMP-quark and WIMP-gluon interaction operators. 
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In the SI part, contributions from individual nucleons in the nucleus add coherently 
and the finite size effects are accounted for by including the SI nuclear form factor F(q). 
The differential cross section for the scalar part then takes the form p[ 

^ = — [Zf p + (A — Z)f n } 2 F\q), (3.2) 

where f p and f n are the effective neutralino couplings to protons and neutrons, respectively. 
Explicit expressions for the case of the supersymmetric neutralino can be found, e.g., in p8|| . 
The formalism we follow has been reviewed in several recent papers [||, |28|, We have 
re-done the original complete calculation of Drees and Nojiri and agreed with their 
results. 

A convenient quantity which is customarily used in comparing theory and experimental 
results for SI interactions is the cross section a p for WIMP elastic scattering of free proton 
in the limit of zero momentum transfer: 

4 1 = ^Ifp ( 3 - 3 ) 

where [i p is defined similarly to fiA above. The analogous quantity for a target with nuclei 
with mass number A can then be expressed in terms of a p T as 

°f = 'A [Zf P + (A- Z)f n f = (^) 2 A*a s /. (3.4) 

One can do so because, for Majorana WIMPs, f p ~ f n . 
The coefficients f p ^ n can be expressed as |L0| 

m v m Q 27 m„ 

* q=u,d,s c,b,t 

where = 1 - J2 q =u,d,s fr] > and frq is S iven b y < P\ m gQl\P >= m P fTq (.9 = 
and analogously for the neutron. The masses and ratios B q = (p\qq\p) of light constituent 
quarks in a nucleon come with some uncertainties. For definiteness, we follow a recent 
re-evaluation [12] and assume m u /md = 0.553 ±0.043, m s /m^ = 18.9 ±0.8, and B(i/B u = 



0.73 ± 0.02, as well as 

/W = o.020 ± 0.004, f^ d = 0.026 ± 0.005, = 0.118 ± 0.062 

fW> = 0.014 ± 0.003, f^} = 0.036 ± 0.008, f!$ = 0.118 ± 0.062, 

which numerically gives SI cross section values very similar to using the set of [Q]. Some 
other recent studies use a new determination of ct^n to derive a much larger value for 
~ 0.37 [ 29 1 . We find such values somewhat questionable since they imply that the 



strange quark component of the nucleon would be larger than the up and down ones. We 
have numerically checked that using the set of input parameters of ]29|] gives typically SI 
cross section values a factor of six higher than in our case. 
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It is worth noting that, despite several different diagrams and complicated expressions, 
it is the exchange of the heavy scalar Higgs that in most cases comes out to be numeri- 
cally dominant. It is further enhanced when the neutralino is a mixed gaugino-higgsino 
state p0[ . As mjj ~ tua increases, a^ 1 drops as m^ 4 because of the t-channel propagator 
effect in XQ —> XI elastic scattering. Eventually, at smaller a^ 1 , squark exchange becomes 
important, and even dominant, instead. 



4. Details of the Scan 

As outlined above, we use seven parameters tan/3, M2, /z, tua, fhq, fhr and A t = Ab to 
conduct a careful scan of the general MSSM parameter space. For their allowed ranges we 
take: 

50 GeV < M 2 < 2 TeV 
50 GeV < \n\ < 2 TeV (4 TeV) 
50 GeV < mj < 2 TeV (4 TeV) 
200 GeV < fhq < 2 TeV (4 TeV) (4.1) 
90 GeV < m A < 2 TeV 

< \At, h \ < ITeV 

5 < tan/3 < 65 

while we set A T = since we treat the masses of the slepton as independent parameters 
anyway. In addition to a general scan of the parameter space, in many cases we do several 
focused scans and explore the effect of extremely large values of fi, fhr and fhq beyond 
2 TeV (given in brackets above) in order to derive parameter-independent lower limits on 
(Tp for a big range of large m x , as described below. The minimum values are set so that 
the resulting physical masses of Higgs and superpartners are limited from below by collider 



bounds. For the lighter chargino we take m ± > 104 GeV [31], for sleptons the lower 



limit of 90 GeV [31| and for squarks 200 GeV [32]. The lower limit on m x depends not 
only on a model but also on a number of additional assumptions |33|. For this reason, 
in our analysis we conservatively do not impose a direct experimental limit on m x , but 
instead infer it from the other limits, especially the one on the chargino mass. As regards 
the lightest Higgs mass m^, in much of the parameter space (jua > 120 GeV) the lower 



limit on the Standard Model Higgs of 114.1 GeV [34j applies. However there are two 
important points to note. Firstly, theoretical uncertainties in computing in the MSSM 
are estimated at 2 — 3 GeV. Conservatively, we thus require only > 111 GeV. Secondly, 
for 90 GeV < rriA < 120 GeV, there still remains a sizable range of the (m/j, m^-plane 
where the lightest Higgs mass given roughly by > 0.78(m,A + 21.7 GeV) is allowed 



We will comment below on the effect of this low-mass range on increasing the largest 
allowed values of . 

Among indirect limits on SUSY, b — > sj often places an important additional constraint 
on the allowed parameter space. We calculate the SM contribution to BR(i? — > X s j) at 
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the full NLO level and include dominant tan /3-enhanced NLO SUSY and also include 
the c-quark mass effect on the SM value @. At this level of accuracy the SM prediction 



is BR(B — > X s j) = (3.70 ± 0.30) x 10 4 [37|. This range is partially overlapping with the 
new world-average |}9| 

BR(B -> X s7 ) = (3.41 ± 0.36) x 10" 4 (4.2) 

which has gone up from the previous range of (3.23 ±0.72) x 10" 4 following the new result 
from BaBAR As described in more detail in [111], with an update in p8j , we accordingly 
allow the full SM+SUSY contribution to be in the range BR(5 -> X s ^) = (3.41 ± 0.67) x 
10 -4 . It is important, however, to stress here an important salient point. In computing the 
SUSY contribution to BR(i? — > X s ^f) one usually makes an implicit assumption of minimal 
flavor violation in the down-type squark sector, which is theoretically poorly justified. 
Even a slight modification of the assumption often leads to a significant relaxation of the 
bound from b — > s'j to the point of even allowing [i < and relatively light superpartner 



masses [39] 



A very recent measurement of the anomalous magnetic moment of the muon = 
(g^ — 2)/2 [jl9| has confirmed a previous value [40 but with twice-increased precision [19] 



a^ xpt - 11659000 x 10~ 10 = (203 ± 8) x 10~ 10 . (4.3) 
The LO hadronic vacuum polarization contribution has recently been re-evaluated in EL 



42] and the light-by-light corrections in ]43| . In [41] the updated SM prediction of [44] has 
been found to be af* - 11659000 x 10~ 10 = (169.1 ± 7.8) x 10 -10 when applying data from 
e + e~ annihilation cross sections and (186.3 =L 7.1) x 1CU 10 when applying r-decay data. 
This leads to a 3 a discrepancy 

Aa M = a^ xpt - aJ M = (33.9 ± 11.2) x 10" 10 (4.4) 

when using the e + e~-based data, or to a 1.6 a deviation 

Aa M = a^ xpt - aJ M = (16.7 ± 10.7) x lO" 10 (4.5) 

when applying the r-based data. 

If interpretted in terms of SUSY, eq. (4.4) restricts the allowed SUSY contribution to 



22.7 x 10~ 10 < a^ USY < 45.1 x 10~ 10 (la) (4.6) 
11.5 x 10~ 10 < a^ USY < 56.3 x 10~ 10 (2a). (4.7) 



Similar ranges are obtained by using the results of [^, 45 1. The eqs. (4.4) and ( f4.5| ) further 



imply that /i > (the sign of the SUSY contribution is the same as that of fj). 

This is clearly an intriguing hint for "new physics". However, since some other recent 
evaluations tend to give a larger value of a^ M |46[[ and a larger error bar [47], at this point 



we will not strictly impose the (g — 2)^ constraint on the parameter space that is otherwise 
allowed by all other constraints. Nevertheless, below we will discuss the important impact 
it has on the upper bounds on m x . 
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Figure 1: Ranges of a^ 1 in the general MSSM vs. m x for tan/3 = 35, A t = Af, = 1 TeV and 
u > 0, which are allowed by the bounds from colliders, b — ► and 0.1 < fl x h 2 < 0.2, but not 
from (g — 2) M . Also marked are some results of recent experimental WIMP searches. The thick 
black line (and a left-pointing arrow) indicates a parameter-independent lower bound on 1 for 
550 GeV < m x < 1020 GeV. No similar bound can be set for lower m x because of the neutralino- 
slepton coannihilation effect, as explained in the text. 



As regards the WIMP relic abundance, a lower limit on the age of the Universe con- 
he* 



servatively gives Q x h 2 < 0.3, while "direct" measurements of the CDM lead to 0.1 < 



n x h 2 < 0.2 which we will treat as a preferred range. Recent studies of the CMBR seem 
to imply even more restrictive ranges; for example, Q x h 2 < 0.12 ± 0.04 in pSj ], We will 
not apply this narrower range yet but will comment below on its impact on the upper 
and/or lower limits on . We compute Q x h? as accurately as one reasonably can, at the 
level of a few per cent both near and further away from poles and thresholds, by applying 
our recently derived exact analytic expressions for neutralino pair-annihilation jllj] and 
neutralino-slepton coannihilation and by using an exact procedure for the neutralino 
coannihilation with chargino and next-to-lightest neutralino 51, 52]. 



5. Results 



The allowed ranges of the SI cross section that result from our scans are illustrated in 
Fig. H for tan/? = 35, A t = = ITeV and fx > 0. Since dp 7 generally grows with tan/3 
due to an enhancement in the heavy scalar Higgs coupling to down-type quarks, the above 
choice is a reasonable compromise between the low and the very large values of tan (3. In 
deriving the allowed ranges of a^j 1 we have imposed the bounds from colliders and from 
b — > sj, and have further required 0.1 < £l x h 2 < 0.2, as described earlier, but not yet the 
bound from (g — 2)^. We also mark with a thick solid line a parameter-independent lower 
bound on cr^ which will be explained below. 

We can see a big spread of a^ 1 over some seven (four) orders of magnitude at small 
(large) m x . In fact, the upper values of a^ 1 exceed the latest experimental limits, including 
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Figure 2: Sensitivity of the upper and lower limits on 1 in Fig. |l| to various assumptions and 
constraints. Upper left window: the light blue region would be allowed if Q, x h 2 < 0.1. The dark-red 
region would be excluded if one neglected the effect of neutralino-slepton coannihilation. Upper 
right window: the light blue region would be allowed if one lifted the constraint from b — > 57. 
The regions to the right of the vertical dotted (dashed) lines are excluded by imposing current 1 a 
(2 a) CL bound from (g — 2)„. Lower left window: the upper dark-red region would be excluded 
by assuming m./, > 111 GeV for all rriA (i.e., by neglecting a window of lighter to/, which is still 
allowed for tua < 120 GeV). Also shown in this window is the effect of restricting itia < 1 TeV. 
Lower right window: the same as for vtla but for /j. 



the recent result from Edelweiss [|| and the new limit from the UKDMC Zeplin I detec- 
tor |5J. Also shown is the CDMS bound and the so-called DAMA region. It is clear that 
today's experiments have already started probing the most favored ranges of cr^ that come 
from SUSY predictions for neutralino cold dark matter. 
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In the four windows of Fig. ^ we show the effect of the most important constraints on 
the upper and lower limits on the allowed ranges of a~ . Firstly, in the upper left window 
we show the effect of relaxing the cosmological bound by allowing Q x h 2 < 0.1. Obviously, 
larger ranges of a^ 1 now become allowed since an enhancement in the neutralino pair- 
annihilation cross section often, by crossing symmetry, implies an increase in (J- 1 ■ Note, 
however, that a combination of all the other constraints, most notably a lower limit on 
rrih from LEP and the constraint from b — > sj, prevents from rising by more than 
about one order of magnitude and only for not very large values of m x . On the other hand, 
imposing a narrower range 0.08 < Q x h 2 < 0.16 has almost no effect on the upper and lower 
limits on a^ 1 , although it does remove a number of points from the allowed ranges of <Jp ■ 

In the same window we also show an important effect, already pointed out in pl| , 
of including neutralino coannihilation with sleptons (predominantly with the lighter stau) 
on allowing very low ranges of a^ 1 at smaller m x . At lowest m x < 120 GeV the relic 
abundance can be reduced to the favored range by choosing in the scan light enough 
sleptons, even without coannihilation. By simultaneously choosing large enough heavy 
Higgs and squark masses, one can reduce to very low values of a few xl0 _12 pb. As 
m x increases, Q x h 2 would normally increase as well, and become too large, but it is there 
that neutralino-slepton coannihilation kicks in. Since dp 1 is independent of the slepton 
masses, by carefully scanning the parameter space, one can always find fhj not much above 
m x , in which case Q x h 2 can be sufficiently reduced again to fall into the favored range. 
The effect is very strong for smaller m x , thus explaining a sharp rise of the left side of the 
dark-red region allowed by neutralino-slepton coannihilation, but, as the process becomes 
increasingly inefficient at larger m x [27, [5(|, it gradually fades away. 

In the upper right window of Fig. ^ we present the effect of imposing the constraint 
from (g — 2)^. We can see that, for this case, m x < 245 GeV (la CL) and m x < 
420 GeV (2 a CL). This upper limit comes from the fact that, as m x increases, the SUSY 
contribution from the x ~ A* an d X~ ~ loops become suppressed and at some point be- 
comes too small to explain the apparent discrepancy between the SM and the experimental 
measurement |5^, [54|]. On the other hand, the upper and lower limits on a^ 1 are not really 
affected. 

In the same window we also show the effect of relaxing the constraints from b — > 57. We 
can see that if it were not imposed, at large m x the upper limit on a^ would significantly 
increase. In this region, the mass of the pseudoscalar Higgs, and therefore also the heavier 
scalar, is rather small, thus giving larger a^ . However, the mass of the charged Higgs is 
then also on the lower side, and a cancellation between a charged Higgs-top quark loop 
and chargino-stop loop contribution is not sufficient to reduce BR(5 — > X s j) to agree with 
the experimental limit. 

However, we remind the reader that a slight relaxation of the underlying assumption 
of minimal flavor violation in the squark sector often leads to a significant weakening of 
the bound from b — > sj to the point of even allowing fi < [39|. We would therefore be 
cautious in applying the b — > 57 constraint rigidly. 

In the lower left window of Fig. [2| we present the sensitivity of the upper limit on a^ 1 
to the lower limit on the light Higgs mass. We can see that it would sizably decrease 
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Figure 3: Sensitivity of to tua and /i for the case of Fig. [l]. We concentrate on the region 
of parameter space where m x ~ 800 GeV. The whole marked region is consistent with all the 
constraints from colliders and b — > S7. By further imposing the constraint 0.1 < i7 x /i 2 < 0.2 one 
selects only the red region. Near the resonance tua — 2m x , significantly smaller values of cr^ 
become allowed by 0.1 < fl x h 2 < 0.2, (mostly in the increasingly pure bino region), but become 
eventually limited from below, independently of increasing the maximum allowed value of fj,. 

if we neglected the region of small 90 GeV < < 111 GeV which is still allowed at 
niA < 120 GeV, and instead required mh > 111 GeV for all tua- Note also that the new 
experimental limits on are for the most part inconsistent with the possibility of the 
light Higgs scalar. 

In the same window and in the lower right window we explore the existence of the lower 
limit on a^ 1 and its dependence on the assumed upper limit on tua and [a, respectively. As 
we can see, the lowest values of a^j 1 are often to a large extent determined by a somewhat 
subjective restrictions from above on these parameters. As one allows either fx or tua above 
1 TeV the lower limit on a~ relaxes considerably. 

However, we argue that, by requiring sizable enough Sl x h 2 {e.g., £l x h 2 > 0.1), it is 
possible to set a parameter-independent lower bound on o~p for a considerable range of large 
m x ~ 1020 GeV (marked with a thick solid line and left-pointing arrow in Fig. ^). The limit 
holds for an arbitrary case of the neutralino (gaugino, higgsino or mixed), independently of 
how large /i and other SUSY parameters are taken. Let us first consider the gaugino limit 
/i 3> M2. In this case the origin of the parameter-independent bound is displayed in Fig. || 
where we plot 1 as a function of the pseudoscalar Higgs tua for m x ~ 800 GeV and scan 
over all the other parameters. As vtla and other parameters are varied, large ranges of a^ 1 
remain allowed by collider and indirect constraints, depending on the maximum allowed 
value of [i. For each fixed [m, o-p 1 decreases proportionally to the fourth power of vtlh — rriA 
because of the (typically dominant) i-channel exchange of the heavy Higgs, as the marked 
cases of \i in Fig. ||| clearly demonstrate. As /i increases, at fixed m x one moves deeper into 
the gaugino region and typically finds large £l x h 2 > 0.2. By imposing 0.1 < Q x h 2 < 0.2 
one selects only a narrow red (dark) range with a large bino component. The bino purity 
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(pg = N^) increases with increasing fx, which normally quickly gives too large Q x h 2 . 
However, for each m x one can choose uia — 2m x (roughly tha = 1600 GeV in Fig. ||) in 
which case Q x h 2 is reduced to an allowed level by a wide resonance due to A-exchange. This 
leads to allowing a much reduced , while still being consistent with 0.1 < £l x h 2 < 0.2. 
However, because of the finite width of the ^4-resonance, at large enough \x one reaches the 
lowest value of <jp (in this case 2 x 10 -11 pb) which is parameter-independent. Away from 
the resonance (for example, if one imposed itia < 1 TeV), one would obtain the lower bound 
Up 1 > 3 x 10 -8 pb, basically independently of whether /i < 1 TeV is imposed or not. In 
fact, one can see this effect in the lower left window of Fig. § where imposing vtla < 1 TeV 
causes the lower limit on a'p 1 to suddenly jump up at around m x ~ m^/2 — 500 GeV. 

The above parameter-independent lower limit arises in the gaugino case and applies 
to arbitrarily large m x . The case of the gaugino-like LSP can be argued to be more 
attractive as being less fine-tuned than the higgsino-like one [26|. In the MSSM the LSP 
is a nearly-pure higgsino for M<i > 300 GeV which, by applying the assumption of gaugino 



mass-unification (see below ( |2.1[ )) to the gluino mass, implies rrig = §^-^2 ^ ITeV and 
therefore large and less "natural" soft SUSY-breaking scale. Nevertheless, in the spirit of 
generality, we need to extend the analysis to the case of the higgsino-like neutralino. 

For the higgsino-like LSP (M2 [i — rn x ), if we remain within the ranges of param- 
eters given in (|4.1| ), we find o~p some two orders of magnitude larger than in the gaugino 
case presented above. However, one can in principle reduce o~p 1 to arbitrarily small values 
by going to the limit of pure enough higgsino {M2 in the multi- TeV range) , in which case 
the neutralino-Higgs coupling would be arbitrarily reduced, and by further suppressing the 
squark contribution by making them extremely heavy. It is therefore reasonable to ques- 
tion the existence of the lower bound on Gp . However, in the multi-hundred GeV range 
(m t < m x < 1020 GeV) Q x h 2 remains typically very small Q x h 2 <C 0.1 although it does 
increase with m x . For tan/3 = 35 it reaches 0.1 for m x < 1020 GeV, 0.2 for m x < 1.6 TeV 
and 0.3 for m x < 2.5 TeV. (These values decrease somewhat with increasing tan/3.) Since, 
as mentioned above, we impose 0.1 < £l x h 2 (or a similar sizable lower limit on Q x h 2 ), in 
the range of m x < 1 TeV displayed in Figs. [I] and ^ the parameter-independent lower limit 
on (jp 1 holds. (In order to be clear that the parameter-independent lower bound on o~p 
applies to a general neutralino so long as m x < 1020 GeV, in Figs. [I] and ^ we have put a 
left-pointing arrow at m x = 1 TeV.) 

For the higgsino-like LSP in the mass range 1020 GeV < m x < 1600 GeV the relic 
abundance is 0.1 < £l x h 2 < 0.2, basically indpendently of how large M2 is, and no lower 
bound on cjp can in principle be set. The range of m x increases to 2.5 TeV if we allow 
& x h 2 < 0.3. Larger values of m x for an arbitrary neutralino gaugino/higgsino composition 
are inconsistent with £l x h 2 < 0.3. This upper limit on m x relaxes to 2.6 TeV for tan (3 = 10. 

On the other side, at low enough m x coannihilation with sleptons prevents one from 
deriving a firm lower limit on cr^ . Indeed, by suitably choosing the slepton mass not too 
much above m x , we can always reduce £l x h 2 below 0.2. Thus, it is possible to set firm lower 
limits on but only for large enough m x and even though they correspond to extremely 
large values of fi and accordingly involve much fine-tunning. 
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Figure 4: The same as in Fig. |] but for tan/? = 10 (left window) and 50 (right window). Also 
marked is the effect of imposing /i < 1 TeV. 



The above discussion of the parameter-independent lower limit on cr^ and upper limit 
on m x has been presented in the case of fi > but it obviously applies also to the case 
fi < 0. 

The dependence of on At and A}, is rather weak. As the tri-linear terms devi- 
ate from zero, the mass of the lightest Higgs generally increases due to somewhat larger 
mass splittings among the stops and sbottoms. As a result, the normally subdominant 
contribution to from the i-channel /i-exchange is slightly reduced. 

A far more important effect is that the number of SUSY configurations satisfying all 
experimental constraints, especially that from b — > sj, decreases significantly. This is 
because the cancellation between charged Higgs loop and chargino loop contribution to the 
b — ► S7 decay rate becomes more inefficient as At decreases. For example, for At = — 1 TeV 
only a handful of points remain allowed. Generally, we have concluded that the regions 
allowed by lower values of At and At, fall into the regions allowed by the choice At = 1 TeV. 

In order to display the dependence on tan /3, in the left and right window of Fig. || we 
present the cases of tan/3 = 10 and 50, respectively. Note that, for small tan/3 = 10 the 
largest allowed values of a^ 1 are roughly an order of magnitude smaller than at tan (3 = 50 
because of the tan /3-dependence of the heavy scalar coupling to down-type quarks, as 
mentioned earlier. Notice a significant decrease in the upper ranges of at large m x for 
tan/3 = 10, which is caused by exceeding the upper limit (4.08 x 10~ 4 ) of the allowed range 
of BR(£? — ► X s j). It is clear that the constraint is more severe in the case tan/3 = 10 rather 
than at larger tan (3. This may sound somewhat counter-intuitive since, for example, in the 
Constrained MSSM, the constraint from b — > sj on the (CMSSM) parameter space becomes 
more pronounced at larger tan (5. This is because, in the CMSSM the pseudoscalar Higgs 
mass, hence also the charged Higgs mass, is typically large, and the corresponding charged 
Higgs-top quark loop contribution to BR(B — > X s ^) becomes suppressed. At smaller 
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Ranges of elf 7 in the general MSSM 



for /i > 0, which are allowed by collider 



S7 and 0.1 < il x h 2 < 0.2. Also marked are some results of recent experimental WIMP 



searches. The thick black line and a left-pointing arrow indicate a parameter-independent lower 
bound. The region below the dashed line is excluded if one imposes the constraint /j < 1 TeV. The 
ranges of m x to the vertical lines are excluded at 1 a and 2 a CL by the current discrepancy between 
the experimental value of (g — 2) M and the Standard Model prediction. 



m\j2 (thus also m x and m x ±) and tuq (thus also light enough stop) the tan /^-dependent 
(negative) contribution from the chargino-stop loop gives too small BR(i? — > X s j), below 
the lower experimental limit, thus producing a strong lower bound on m x at not too large 
mo- I n contrast, in the general MSSM case, we can choose small values of tjia which is 
a free parameter. This small txia implies a big positive charged Higgs loop contribution 
to BR(B — ► X s j). For smaller tan/3 and at smaller m x this is reduced to an acceptable 
range by the chargino-stop loop contribution. However, at large m x , the chargino-stop 
loop cannot cancel the charged Higgs contribution anymore and one exceeds the upper 
experimental limit on BR(i3 — > X s ^). Since the chargino loop contribution is proportional 
to tan/3, at large tan/? the cancellation can be achieved even at large m x . For this reason, 
in the right window of Fig. || there is no analogous decrease in the largest allowed cr^ at 
large m x , in contrast to the left window. Clearly, the constraint from b — > 57 is more severe 
for smaller tan (3. 

A similar effect of can be observed in the case of the — constraint. By comparing 
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Fig. ^| with the upper right window of Fig. ^ one can see that it produces as stronger upper 
bound on m x at smaller tan (3. For example, imposing a la bound implies m x < 120 GeV 
for tan/3 = 10, while for tan/3 = 50 the bound moves up to m x < 310 GeV, as denoted in 
the two windows of Fig. Also marked is the effect of imposing \x < 1 TeV. While not 
being a firm constraint, it does, in our opinion, indicate the region which may be considered 
as somewhat less fine-tuned. 

As discussed above, for a considerable range of large m x < 1020 GeV (marked with 
a thick solid line and left-pointing arrow in Fig. ^), we can again set up the parameter- 
independent lower limit on <7p , in analogy with the case tan (3 = 35 of Fig. [|. In order to do 
so, we had to explore extremely large ranges of /i up to some 4 TeV at smaller tan (3 in order 
to saturate the bound Q x h 2 < 0.2. On the other hand, at smaller m x < 440(800) GeV 
for tan/3 = 10(50) the lower limit on remains basically independent of tan/3 since 
£l x h 2 at lower m x is determined mostly by the coannihilation with sleptons, as discussed 
in detail in the case tan (3 = 35. In this case the lower limit on a^ 1 arises from restricting 
H below the value for which the parameter-independent lower limit arises at larger m x . 
Notice that the thick line extends to lower m x at smaller tan/3 because the efficiency of 
the neutralino-slepton coannihilation increases with tan (3. 

Finally, in Fig. || we summarize the results for the full scan conducted so far for tan (3 = 
5, 10, 35, 50, 55, 60, 65 and for fj, > 0. We repeat that, in determining the allowed (blue) 
region we applied the constraints from collider searches, b — > 57 and 0.1 < {l x h 2 < 0.2. 
The lower limit on mostly comes from low tan (3 and, where applicable, we mark with 
the solid line where it is parameter-independent. The effect of restricting /i < ITeV is 
marked with a dashed line. Also indicated is the impact of the new measurement of (g — 2)^ 
on the mass of the neutralino. If confirmed, the e + e~-based range ( [4.4| ) will imply rather 
stringent upper limits on m x 

m x < 350GeV (IctCL) (5.1) 
m x <510GeV (2 a CL). (5.2) 

If the r-based numbers (t4.5|) are applied instead, one obtains m x < 800 GeV (1 a CL) and 
no upper bound at 2 a CL (for fi > 0). 

The vast ranges of a^ 1 predicted in the framework of the general MSSM may be 
somewhat discouraging to DM WIMP hunters. It is worth noting, however, that it is the 
region of smaller m x , below a few hundred GeV, that not only is implied by the new result 
for (g — 2)^, but is also theoretically more favored as corresponding to less fine-tuning. 
Furthermore, ranges of very small 10~ 12 pb < <jp < 10~ 8 pb generally correspond either to 
very large (and therefore perhaps somewhat less natural) values of fi and/or m^, or become 
allowed by selecting slepton masses on the light side, and in the x~sl e Pton coannihilation 
region, within some 20 GeV of m x , which again can be be considered as a finely-tuned 
case. 

At the end, we comment again on the case of the Constrained MSSM. Because the 
model is much more restrictive, the ranges of cr^ that one obtains in the parameter space 
allowed by all constraints, are very much narrower 12, |l4|, [l5|, |l6[| . They are also 
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typically somewhat lower than the largest ones allowed in the general MSSM. For example, 
at tan/3 = 50 we find a^ 1 ~ l(T 7 pb at m x = 100 GeV and a^ 1 ~ 7 x 10~ n pb at the 
largest (neglecting (g — 2)^) allowed value of m x = 800 GeV. On the other hand, because 
the model is defined at the grand-unified, and not electroweak, scale, in the case of large 
tan/3 > 50 and/or large scalar masses, theoretical uncertainties involved in the running of 
parameters are substantial and have much impact on the resulting ranges of both m x and 
. We will explore the case of the Constrained MSSM in an oncoming publication. 



6. Conclusions 

We have delineated the ranges of the SI cross section a^ 1 in the general MSSM, which are 
consistent with current experimental bounds and for which one finds the expected amount 
of dark matter. We have further discussed the dependence of our results on the experi- 
mental constraints and on the underlying theoretical assumptions. While the ranges which 
we have obtaine extend over more than six orders of magnitude, we find it encouraging 
that the experimental sensitivity that has recently been reached, now allows one to explore 
our theoretical predictions for the MSSM. As we have argued above, smaller values of the 
WIMP mass and also larger values of may be considered as more natural, which will 
hopefully be confirmed by a measuring a positive WIMP detection signal in the near future. 
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